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ABSTRACT 

We observed three high-mass star-forming regions in the W3 high-mass star formation complex 
with the Submillimeter Array and IRAM 30 m telescope. These regions, i.e. W3 SMS1 (W3 IRS5), 
SMS2 (W3 IRS4) and SMS3, are in different evolutionary stages and are located within the same 
large-scale environment, which allows us to study rotation and outflows as well as chemical properties 
in an evolutionary sense. While we find multiple mm continuum sources toward all regions, these 
three sub-regions exhibit different dynamical and chemical properties, which indicates that they are 
in different evolutionary stages. Even within each subregion, massive cores of different ages are found, 
e.g. in SMS2, sub-sources from the most evolved UCHii region to potential starless cores exist within 
30 000 AU of each other. Outflows and rotational structures are found in SMS1 and SMS2. Evidence 
for interactions between the molecular cloud and the Hn regions is found in the 13 CO channel maps, 
which may indicate triggered star formation. 

Subject headings: stars: formation - stars: massive - ISM: jets and outflows - ISM: molecules - ISM: 
individual (W3) 



1. INTRODUCTION 

The formation and evolution of high-mass stars is 
still poorly understood. High-mass star formation not 
only proceeds exclusively in a clustered mode, but it 
becomes increasingly apparent that several episodes of 
high- and low-mass star formation can proceed sequen- 
tially within the same molecular cloud comple x. Well- 
know n examples are NGC 63341 fc I(N) (e.g., iSandelll 

J [2005; 



20001) or the S2 55 complex (e.g., iMinier et all 



Wang et al.ll201l| ). Since the general molecular and stel- 
lar environment in these regions is almost the same 
within their different subregions, the environmental ef- 
fects are minimized, which allows us to study the high- 
mass star formation in an evolutionary sense. Al- 
though there has been considerable progress in the gen- 
eral understandin g of high-mass star formation over the 
last d ecade (e.g. JZinnecker fc Yorkdl2007c IBeuther et al.l 
2007), many questions still remain open. The time evo- 
lution of the chemical and physical properties (e.g., ro- 
tating/infalling envelope, outflow) is still unclear. 

Among the physical properties, massive molecular out- 
flows are among the most studied phenomena. The out- 
flows have been detected in various high-mass star forma- 
tion stages, from the early stage as IRDCs to relatively 
evolved high- mass hot cores (e.g., IBeuther et al.l 12002: 
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I Wu et all 120041: IZhang et al.l 120051: IBeuther fc Nissenl 
120081 : iFallscheer et al-ll2QQ9b IVasvunina et al.ll2011l ). Ob- 
servations towards the S255 complex show that out- 
flows in the very young and cold region have low ve- 
locity and are relatively confined around the driving 
source, while the outflows associated with the relatively 
evolved high-mass protostellar object s exhibit large scal e 
collimated high velocity s tructure (jWang et al.l 120 111 ). 
IBeuther fc Shepherd! (|2005l ) proposed an evolutionary se- 
quence for the massive outflows in the high-mass proto- 
stellar phase: as the central object continues to accrete, 
it evolves from a B star to an O star and the outflows be- 
come less collimated under the pressure of the increased 
radiation and stellar wind. 

Massive disks are a key component for high- 
mass stars to overcome radiatio n pressure and form 
through disk accret ion (e.g., iKrumholz et al.ll2007l 120091 : 
iKuiper et al.l l2Q10h . Rotating disks have been de- 
tecte d around inte r mediat e -mass sources (~ 5 — 1 , , 
e.g., IZhang et all Il998ri ICesaroni et~atl I2005L 120071 1 
but for higher mas s objects, only rota ting structures 
known as "to roids" (ICesaroni et al l l2007f) have been de- 
tected (e.g., Beltran et al.l 120041 : IBeuther et "all l2005ri 
IFallscheer et al.l [2009V These structures are not in Ke- 
pleri an rotation and have typical radii of thousands of 
AU. Fallscheer (2010) sees an evolutionary trend that 
the size of the potent ial disk signature decreases from 
IRDCs to hot cores. I Wang etaLl (|2011l ) also detected 
rotating toroids associated with high-mass protostellar 
objects and hot cores in the S255 complex, but the evo- 
lutionary trend is not clear according to their results. 

From IRDCs to high-mass hot cores, massive 
star f ormation exhibits diverse chemi c al properties 
(e.g.. ISutton et all [19851: ISchilke et al.l Il997ai [2 001: 
Beut her et all 120091 : fVasvunin aet al.l 1201 lflWang et al.l 



2011[ ), yet relatively few evolutionary details are known. 
High resolution observations als o show spatial varia- 
tions among many species (e.g., IBeuther et al.l [2005c, 
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12009b iWang et all 1201 Several theoretical studies ex- 
ist on the chemical evolution during massive star forma- 
tion (e.g.. iCharnlev et al.lfT997L IN omura fc Millarl[200l 
IViti et alJl2QQ4l iBeuther et al.ll2009D . however, the high- 
resolution observational database to test these models is 
still relatively poor. 

In this work we present high resolution SMA inter- 
ferometry observations as well as single-dish CO (2 — 1) 
observations of three subregions (W3 SMS1, W3 SMS2 
and W3 SMS3) wi t hin the same environm e nt, W3 Main 
(lLadd et all 119931: [Tieftrunk et al.l 119981 : iMoore et all 
120071 : Megea th et all 12008( 1. We study the evolution of 
the chemical and physical properties (outflows, rotation), 
ranging from potential starless cores to the most evolved 
ultracompact Hn (UCHii) regions. W3 Main is part of 
the W3 molecula r cloud and lo cated at a distance of 
1.95 ± 0.04 koc dXu et all 12006ft . Single dish sub-mm 
observations (Lad d et al.lll993f ) resolved three continuum 
sources, W3 SMS1, SMS2 and SMS3 (Fig.HJ). 

The easern region W3 SMS1: — The central infrared 
source W3 IRS5 is associated w ith the hypercompact 
Hn (HCHu) c luster W3 M (jClaussen et al.l fl99l 
iTieftrunk et al.l 119971 ) with a total luminosity of 
2xl0 5 Lq ([Campbell et al.l 119951 : IWvnn- Williams et al.l 
I1972D. and co in cides w ith water a nd OH masers 
(Clausse neTall 119941: ISarma et al.l I2QQ1L 120021 : 
llmai et al.l 120001 : iGaume fc Mutell I1987D . The low 
radio luminosity of these HCH11 regi ons suggest s they 
might be stellar win d /jet so urces dHoare et al.l 2007; 



iHoare fc Franco! 120071 ). and IWilson et al.l (j2003D found 



two of the HCH11 regions show large proper motion and 
might be jet-knots expelled from the massive stars in 
IRS 5. In the high resolution (< .1") radio observations 
done by Ivan der Tak et all ((2005) these sources do not 
show elongated morphology and the authors argue that 
there are three HCH11 regions together wi th several 
shock-ionized clumps. Meg eath et al.l (|2005f ) reported 
seven near-IR sources with HST observations, and 
three of them have both mid-IR and r adio continuum 
counterparts (|van der Tak et all |2005[). Usi n g the 
Plateau de Bure Interferometer, iRodon etaLl J2008) 
resolved the main continuum source into five sources. 
Bipolar o utflows and even mul ti ple outflows a re also 



reported ( Elaussen et al.l [T9 84: MitcheinFal] 1992 



]hoi et a l. 1993; Hasega wa et al.l 119941 llmai et al.l 120001 : 
IRodon et al.l 2008). It has been proposed that IRS5 is a 
proto- Trapezium and might emerge as a bound Trapez- 
ium similar to that in th e Orion nebula ([M egeath etahl 
2005: IRodon et al.l 120081 ). iMegeath et "all (|1996l ) iden- 
tified an embedded cluster of 80—240 low-mass stars 
centered on IRS5, which in tu rn is surrounded by 
hund reds of low-mass stars (Feigelso n fc Townslevl 
l2008h . W3 IRS7 (associated with the UCHii region W3 
F) and IRS 3 (associated w ith the compact H11 region 
W3 B) are also found here (IWvnn- Wil liams e fall 119721 : 
Harris fc Wvnn-Williamsl 1 19761 : ITieftrunk et alTll997l ). 

The northwestern region W3 SMS2: — The infrared source 
W3 IRS4 lies at the center of the continuum sourc e 
with a luminosity of 6xl0 4 L© (jCampbell et al.l fl995). 
IRS4 is classified as an O8V-B0.5V star and is thought 
to be the excit i ng source of the UCH ii region W3 C 
(jBik et al.ll20lH : iTieftrunkeT al. 1997). The HCH11 re- 



gion W3 Ca lies east of IRS4 which indicates that a mas- 
sive star is forming. OH masers and evidence for an 
outflow are also found tow ards IRS4 (jGaume fc Mutel 
119871 : lHasegawa et al.l[l994fl . 

The southwestern region W3 SMS 3: — W3 SMS3, located 
south of W3 IRS4, is consider ed a quiescent region in 
W3 Main (jMegeath et al.ll2008D . It exhibits strong mm 
continuum emission (Fig. [T]) but no other signs of active 
star formation such as n ear infrared or mid infr ared emis- 
sion (|Ruch et al.ll2007D . ITieftrunk et~aTl (|1998f ) found ex- 
tended emission towards SMS3 with VLA NH3 observa- 
tions, which show quite narrow line widths towards this 
region. ITieftrunk et al.l (pL998) also found four compact 
NH3 clumps with masses around 30 M surrounding the 
main clump. 

Furthermore, NIR imaging and spectroscopy studies 
by iBik et al.l (|2011f ) revealed an age spread of at least 
2-3 Myr in W3 Main and suggested a sequential star 
formation process. While W3 SMS1 and SMS2 show 
signs of active star formation and are surrounded by a 
NIR cluster, SMS3 is still quiescent with no signs of on- 
going star formation. Therefore, W3 Main is an ideal 
region to simultaneously investigate several sites of mas- 
sive star formation at different evolutionary stages within 
the same over all environment. 

2. OBSERVATIONS AND DATA REDUCTION 

2.1. Submillimeter Array observations 

The W3 Main complex was observed in three fields 
with the Submillimeter Array 8 (SMA) on December 31st 

2009 in the compact configuration and on January 30th 

2010 in the extended configuration, both with seven an- 
tennas. The phase centers of the observations, which 
are known as W3-SMS1, W3-SMS2, and W3-SMS3, were 
R.A. 02 h 25 m 40.68 s Dec. +62°05'51.5" (J2000), R.A. 
02 h 25 m 31.22 s Dec. +62°06'25.5" (J2000), and R.A. 
02 h 25 m 29.49 s Dec. -f 62°06'00.6" (J2000), respectively. 
The SMA has two spectral sidebands, both 4 GHz wide 
and separated by 10 GHz. The receivers were tuned to 
230.538 GHz in the upper sideband (v\ sr = —43 km s _1 ) 
with a spectral resolution of 1.2 km s _1 . For the com- 
pact configuration, the bandpass calibration was derived 
from the observations of the quasar 0854+201. Phase 
and amplitude were calibrated with regularly interleaved 
observations of the quasar 0136+478 (15.8° from the 
source). Since the flux measurements of the gain cali- 
brator 0136+478 on December 31st 2009 are available in 
the SMA calibrator database, we do the flux calibration 
with this quasar, and the flux scale is estimated to be ac- 
curate within 20%. For the extended configuration, the 
bandpass calibration was derived from the quasar 3C273. 
Phase and amplitude were calibrated with regularly in- 
terleaved observations of the quasar 0102+584 (11° from 
the source). The flux calibration was derived from ob- 
servations of Titan, and the flux scale is estimated to 
be accurate within 20%. We merged the two configu- 
ration datasets, applied different robust parameters for 
the continuum and line data, and obtained synthesized 

8 The Submillimeter Array is a joint project between the Smith- 
sonian Astrophysical Observatory and the Academia Sinica Insti- 
tute of Astronomy and Astrophysics and is funded by the Smith- 
sonian Institution and the Academia Sinica. 
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beam sizes between 1.5" x 1.1" (PA -35°) and 2.2" x 1.8" 
(PA -34°). The 3a rms of the 1.3 mm (225 GHz) con- 
tinuum image is ~ 10.8 mJy beam -1 and the 3a rms 
of the line data is 0.11 Jy beam -1 at 2 km s -1 spectral 
resolution. The flagging and calibration wer e done with 
the IDL superset MIR (jScoville et al.l fl993l ) which was 
originally developed for the Owens Valley Radio Obser- 
vatory and adapted for the SMA 9 . The imaging and dat a 
analysis were conducted in MIRIAD ([Sauft et al.| [l995). 

2.2. Short spacing from the IRAM 30 m 

To complement the CO (2 — 1) observations which 
lack the short-spacing information and to investigate the 
large-scale general outflow properties, we observed the 
three subregions with the HERA array at the IRAM 30 m 
telescope on February 21st, 2011. The 12 CO(2 - 1) line 
at 230.538 GHz and the 13 CO(2 - 1) at 220.399 GHz 
were observed in the on-the-fly mode. We mapped the 
whole region with a map of 4' x 4' centered at R.A. 
02 h 25 m 41.01 s Dec. +62°07'02.1" (J2000). The sampling 
interval was 0.8". The region was scanned once along the 
declination and right ascension directions, respectively, 
in order to reduce scanning effects. The spectra were 
calibrated with CLASS, which is part of the GILDAS 
software package. The 12 CO data have a beam size of 
11", and the rms noise level of the corrected T m b scale 
is around 0.25 K at 1.2 km s -1 spectral resolution. The 
13 CO data have the same beam size, but the rms noise 
is 0.35 K at 1.2 km s -1 spectral resolution. 

After reducing the 30 m data separately, single-dish 
12 CO and 13 CO data were converted to visibilities and 
then combined with the SMA data using the MIRIAD 
package task UVMODEL. The synthesized beam of the 
combined data is 3.0" x 2.5" (PA -34°). 

3. RESULTS 

3.1. Millimeter continuum emission 

Figure [2] shows the 1.3 mm continuum maps of the 
three regions. We resolve several continuum sources in 
all three regions. Assuming optically thin dust emis- 
sion, we estimate the gas mass and column density of 
th e continuum sourc es fo llowing the equations outlined 
in lHildebrandl (|1983f ) and lBeuther et al.) (|2005cft . We as- 
sume a dust temperature of 40 K, grain size of 0.1 /im, 
grain mass density of 3 g cm -3 , gas-to-dust ratio of 100, 
and a grain emissivity index of 2 (corresponding to n ^ 
0.3 c m 2 g -1 for comparison with Ossenkopf fc Henningl 
1994). Here we assume the same dust temperature for 
all the continuum sources to simplify the calculation. 
The fact that the continuum sources may have differ- 
ent temperatures (see Sect. 13. 5p would also bring uncer- 
tainty to the mass results. The SMA continuum mea- 
surements listed in Table [1] are not corrected for the 
primary beam, since for SMS1 and SMS2, the primary 
beam correction would only change the flux of our de- 
tected continuum sources by < 2%, and for SMS3 the 
flux would be changed by < 10%. For the SMA contin- 
uum sources which are associated with radio emission, 
the free-free emission is subtracted. We assume optically 
thin free- free emission (spectral index ~ —0.1), and use 

9 The MIR cookbook by Chunhua Qi can be found at http: / /cfa- 
www.harvard.edu/~cqi/mircook.html. 



the VLA 22.5 GHz flux (Ti eftrunk et~all ll997) to esti- 
mate the free- free flux at 225 GHz. In the following, we 
label the mm sources according to a scheme consisting of 
the name of the large-scale region (SMS1 to SMS3) fol- 
lowed by the mm peaks within them, e.g., SMS1-MM1. 
The corrected results are shown in Table [TJ For sources 
SMS1-MM3 and SMS2-MM4, the free-free flux we obtain 
at 225 GHz is even larger than our SMA 1.3 mm mea- 
surements, which is likely due to the fact that the VLA 
observations have better uv coverage and suffer less from 
missing flux than the SMA observations. 

With the same set-up, we also calculate the gas mass 
and the column density of the SCUBA 850 jira contin- 
uum peaks shown in Table [2] (The SCUBA data were 
downloaded from the archive, Di Fra ncesco et al.l [2008) 
and compare the results with the ones in Table [TJ The 
masses we obtain from the SMA observations for the W3 
SMS1, W3 SMS2 and W3 SMS3 are only 9%, 12% and 
3% of the SCUBA 850 /xm measurements, respectively. 
For the W3 SMS1 region, we resolve part of the Hn region 
W3 B from the SMA observations (Fig. [2]). However, it 
is already near the edge of our SMA field of view, and 
the emission we recover should be dominated by free-free 
emission. Therefore, we do not include the emission at 
W3 B from the SMA observations for the comparison. 
Since our interferometer observations are not sensitive 
to spatial scales > 27" (corresponding to the shortest 
baseline = 7.5 kA), we filter out the smoothly distributed 
large-scale halo and observe only the compact cores. Fur- 
thermore, the observational fact that we filter out more 
flux for the youngest region W3 SMS3 indicates that at 
the earliest evolutionary stages the gas is more smoothly 
distributed than during later stages where the collapse 
produces more centrally condensed structures. Never- 
theless, all of the column densities we derive are above 
the estim ated threshold for high-ma ss star formation of 
1 g cm -2 (jKrumholz fc Mc Kee 2008), which corresponds 
to a column density of 7Vh 2 ~ 3 x 10 23 cm -2 . 

W3 SMS1— The SMA continuum image of W3 SMS1 
is shown in the bottom left panel of Fig. [2j The 
SMA image reveals four continuum sources, SMS1-MM1, 
SMS1-MM2, SMS1-MM3 and SMS1-MM4 in this re- 
gion. As shown in Table [3j the strongest source SMS1- 
MM1 coincides with the near-infrared s ource W3 IRS5 
and t he HCHh region clu ster W3 M (Clau ssen et al.l 
119941 : iTieftrunk et al.l Il997h . It is also associated with 
the mid-infrared s o urces MIR1 and MIR2 detected by 
van der Tak et all ((2005). H 2 masers are also de- 
tected at the SMS1-MM1 position (jlmai et al.l 120001 : 
Sarmaet al. 2001, 2002). A previous higher resolution 
study shows that SMS1-MM1 is resolved into 4 i ndividual 
1.4 mm continuum sources (|Rod6n et al.l 12008). SMS1- 
MM2 and MM4, which have not been detected before, are 
not associated with any infrared source. S MS1-M M4 i s 
close to an NH3 core detected by Tieft runk et al.l ([1998). 
SMS1-MM3, south of SMS1-MM1, is associated with the 
infrared source W3 IRS7 and the UCH11 region W3 F, 
Table [T] shows that the continuum emission at this region 
is dominated by free- free e mission from the UCH11 region 
and it is relatively evolved (jTieftrunk et al.lfl997h . 

W3 SMS2 — The bottom middle panel of Fig. [2] presents 
the SMA continuum image of W3 SMS2. We detect 7 
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Table 1 

Millimeter continuum peak properties. 



Source 



R.A. 
(J2000) 



Dec. 
(J2000) 



(mjy/beam) 



S v 
(mJy) 



Mass 
(M ) 



free-free l v 



free-free S u 



Nh 2 _ 

(10 24 cm 2 ) measured I u measured S u 



SMS1-MM1 


02:25:40.77 


+62:05:52.3 


123 


405 


40 


6.5 


0.03 


0.02 


SMS1-MM2 


02:25:40.48 


+62:05:50.5 


63 


229 


23 


3.4 


n/a 


n/a 


SMS1-MM3 


02:25:40.54 


+62:05:40.7 


22 


25 




0.3 


0.73 


2.04 


SMS1-MM4 


02:25:39.91 


+62:05:51.7 


22 


35 


4 


1.2 


n/a 


n/a 


SMS2-MM1 


02:25:31.51 


+62:06:22.7 


68 


208 


21 


3.6 


n/a 


n/a 


SMS2-MM2 


02:25:31.36 


+62:06:20.9 


67 


121 


12 


3.5 


0.02 


0.01 


SMS2-MM3 


02:25:31.93 


+62:06:23.9 


55 


113 


11 


2.9 


n/a 


n/a 


SMS2-MM4 


02:25:30.94 


+62:06:15.9 


50 


220 




1.3 


0.52 


1.01 


SMS2-MM5 


02:25:31.05 


+62:06:19.9 


43 


75 


2 


1.3 


0.42 


0.74 


SMS2-MM6 


02:25:31.16 


+62:06:18.1 


34 


50 


0.1 


0.6 


0.69 


0.97 


SMS2-MM7 


02:25:30.85 


+62:06:23.3 


26 


67 


7 


1.4 


n/a 


n/a 


SMS3-MM1 


02:25:28.04 


+62:06:03.6 


21 


40 


4 


1.2 


n/a 


n/a 


SMS3-MM2 


02:25:30.15 


+62:06:01.2 


21 


66 


7 


1.2 


n/a 


n/a 


SMS3-MM3 


02:25:28.21 


+62:06:04.2 


19 


33 


3 


1.1 


n/a 


n/a 



Note. — "n/a" indicates no free-free emission is detected towards this source. For sources SMS1-MM3 and SMS2-MM4, 
the free- free flux we obtain at 225 GHz is even larger than our SMA 1.3 mm measurements, so we do not have the mass 
estimation. 



Table 2 

SCUBA 850 /im continuum data. 



Source 


Iu 


S u 


Mass 






(Jy/beam) 


(Jy) 


(M ) 


(10 23 cm" 2 ) 


W3 SMS1 


19 


39 


715 


5.9 


W3 SMS2 


15 


24 


435 


4.7 


W3 SMS3 


14 


24 


445 


4.4 



Table 3 

Millimetric, radio, NIR and MIR counterparts. 



MM 



radio a 



NIR 



MIR 



SMS1-MM1 W3 M (a+b+c+dl+d2+f) 
SMS1-MM3 W3 F 

SMS2-MM2 W3 Ca 

SMS2-MM4 W3 C 

SMS2-MM5 W3 C 

SMS2-MM6 W3 C 



1+2 C 



IRS5 b 
IRS7 b 

IRS4-b d IRS4-b d 
IRS4-a d 



a fClau ssen et al. (1994); Tieftr unFetaTI (TT997T ) 
b Wvnn- Williams et al. (1972) 
c Ivand er Tak et al. (20051) 
d Figs. [3] and [8] 

continuum sources in this region. As listed in Table [3j 
SMS2-MM2 is associa ted with the HCHu region W3 Ca 
(jTieftrunk et all 11997ft . and Fig. [3] shows that SMS2- 
MM2 also coincides with another elongated infrared 
source IRS4-b which could be tracing the jet driven by 
SMS2-MM2 (see Fig. [13] in Sect. E3J). The elongated 
infrared source IRS4-b is very luminous in Spitzer/IRAC 
bands and dominates the mid-infrared emission here. 
The NI R source IRS4 m arked here as IRS4-a fades out 
(Fig. El iRuch et al.ll2007ft . The top panel of Fig. [2] shows 
that SMS2-MM4, SMS2-MM5, SMS2-MM6 are associ- 
ated with the ultracompact drop-shaped Hn region W3 
C, and Table [1] shows that the continuum emission of 
these three sources is dominated by the free-free emis- 
sion from the UCHii region W3 C. This indicates that 
SMS2-MM4, SMS2-MM5 and SMS2-MM6 are just part 
of the UCHii region W3 C. The other three continuum 
sources, SMS2-MM1, MM3 and MM7, are not associated 



with any infrared or radio continuum source. They are 
likely to be at a very young evolutionary stage. 

W3 SMS3— The SMA continuum image of W3 SMS3 
is shown in the bottom right panel in Fig. [2] The SMA 
observations resolve three continuum peaks in this re- 
gion. These continuum sources are not associated with 
any near infrared or radio source. The Spitzer/MIPS 
24 /im map also shows no emission toward this region 
(|Ruch et al.l 12007). Meanwhile, the gas mass is only 3% 
of the SCUBA 850 /im measurements, which indicates 
that the gas in this region is smoothly distributed in- 
dicative of an early evolutionary stage. All these features 
indicate that W3 SMS3 is an extremely young region. 

3.2. Spectral line emission 

We detect 41 lines from 13 species towards all three re- 
gions. In addition to the three CO isotopologues 12 CO, 
13 CO and C 18 0, we also detect dense/hot gas tracers 
CH 3 CN, CH 3 OH, HNCO, HC 3 N, OCS, shock tracer SiO, 
temperature determinator H2CO, ionized gas tracer hy- 
drogen recombination line H30a, and many SO and SO2 
isotopologues. The detailed chemical properties of each 
region will be explained in the following sections. We 
extract spectra of several molecular transitions at the 
position of each continuum peak within one beam size 
(2.1" x 1.7" for the lower sideband and 2.0" x 1.6" for 
the upper sideband), and the detected lines are listed in 
TableEl Figured shows the spectra towards SMS1-MM1 
and SMS3-MM1 as examples. 

W3 SMS1— In the W3 SMS1 region, we detect 36 
lines from 12 species. Besides the three CO isotopo- 
logues, we also detect many sulfur-bearing species ( 13 CS, 
OCS, SO, SO2 and their isotopologues), some dense gas 
molecules which are usually used t o trace high-mass 
hot cores, such as CH3OH , CH3CN dN omura & Millar 
120041 : iBeuther et all 120091 : iSutton et all 11985ft . and 
the well known kinetic temperature determinator 
H 2 CO {M angum & Wootten 1993; Jansen etld] 119941 : 
iMiihle etalJ 120071 : IWatanabe fc M itchell 2008). All the 
lines we detect have lower energy levels E\ ower /k between 
5.3 to 881 K (Table El). 
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Figure 1. The SCUBA 450 /xm (top) and 850 fim (bottom) con- 
tinuum images. The contour levels are in steps of 5cr starting at 5cr 
in both panels. For the 450 /im panel, a = 0.1 Jy beam -1 , and for 
the 850 /im panel, a = 0.4 Jy beam -1 . The dashed circles mark 
the primary beam of our SMA observations, the filled triangles in 
the top panel mark the SMA continuum sources we resolve, and 
the stars in the bottom panel mark the near infrared sources. 

Figures [5] and [6] present the integrated line maps of all 
species (except SiO, 12 CO and 13 CO which will be dis- 
cussed in detail in Sect. E21). For CH 3 CN(12 fe - ll fe ), 
we only present the line map of k = 0, 1, since the other 
line maps are similar to the one shown. For CH3OH, 
we only present the line maps of two transitions, and for 
H2CO we only show the map for one transition. The 
CH3OH(102,8 — 93,7)A+ line map shows an elongated 
structure closely aligned with the direction of outflow- 
a (see Fig. HI Sect. E3J). All the other CH 3 OH and 
CH3CN lines do not peak exactly on any continuum 
source but somewhere between SMS1-MM1 and SMS1- 
MM2 (Fig. [5j). We derive the peak positions of the 
CH 3 CN(/c = 0,1), CH 3 OH(4 2 , 2 - 3i >2 )E and R30a in- 
tegrated emission through Gaussian fitting, and list the 
results in Table HJ All three H2CO lines show an elon- 



Table 4 

CH3CN, CH3OH and R30a peak positions. 



lines R.A. (J2000) Dec. (J2000) 

CH 3 CN(fc = 1, 2) 2:25:40.54±0.02 +62:05:50.9±0.3 

CH 3 OH(4 2 2 - 3i 2 )E 2:25:40.47±0.01 +62:05:51. 0±0.2 
H30a a ' 2:25:40.68+0.03 +62:05:51.7+0.5 

a We only fitted the peak at SMS1-MM1. 

gated structure which extends in the northwest-southeast 
direction centered between the SMA continuum peaks 
SMS1-MM1 and SMS1-MM2. 13 CS line emission is off- 
set from the SMA continuum peaks (Fig. [5j), similar off- 
set C 34 S emission has also been found bv iBeuther et al.l 
(2009). The hydrogen recombination line H30a is also 
detected towards SMS1-MM1, SMS1-MM3 and the Hn 
region W3 B. Table |4] shows that the H30a peak at 
SMS 1-MM1 is between t he compact radio sou rces a and 
b (jClaussen et al.l 119941 : iTieftrunk etHI 119971 ). and the 
emission might be the joint contribution of sources a, 
b, c, dl/d2 and f (jClaussen et al.l 119941 : ITieftrunk et al.l 
11997( 1. We detect many sulfur-bearing lines in this re- 
gion. Figure [6] shows that all lines of SO and SO2 iso- 
topologues peak on the continuum source SMS1-MM1. 
Furthermore, the 80(65—64) line emission shows an elon- 
gated structure along the red-shifted outflow (outflow-b 
in Fig. [TTJ see Sect. I3.3j) and could be affected by the 
outflow. 

Most of the lines we detected in the SMS1 region are 
found towards SMS1-MM1, except the /c-ladder lines of 
CH3CN and the 13 CS lines. Far fewer lines are detected 
towards SMS1-MM2. Although the hot core tracers 
CH 3 CN and CH3OH are detected toward SMS1-MM2, 
the line emission peaks do not coincide with the contin- 
uum peak of SMS1-MM2. In addition, CH 3 OH emis- 
sion is also associated with the red-shifted outflow of 
SMSl-MMl(outflow-a in Fig. El see Sect. [33]), these 
lines are more likely due to the shock heating excited by 
the outflow. As a typical dense gas molecule, we can- 
not explain why no CH3CN emission is detected towards 
SMS1-MM1. Many lines we label as detected towards 
SMS1-MM2 in Table [6] do not show an emission peak at 
SMS1-MM2 (Figs. |5]and[6j). The reason could be that 
some of the lines associated with SMS1-MM1 are actu- 
ally extended emission and not confined to SMS1-MM1. 
Since we extract the spectra with one beam size, the spec- 
tra we extract towards SMS1-MM2 could be affected by 
the extended emission around SMS1-MM1, which means 
the actual number of lines which originate from SMS1- 
MM2 is even smaller. For SMS1-MM4, only the three 
CO lines, the SO (65 — 5 4) li ne and the SiO line are de- 
tected fFig. lTTj see Sect. 13.31) . SiO is a well-known shock 
tracer (e.g.. IScnilke et al.lll997bf ). and the SiO emission 
detected from SMS1-MM4 is clearly due to the shock ex- 
cited by the outflow. Therefore, SMS1-MM2 and SMS1- 
MM4 appear to be chemically younger than SMS1-MM1. 
No lines are detected towards SMS1-MM3 except CO and 
H30a lines which indicates that this source is relatively 
evolved and has no high density gas in its circumstellar 
environment. 

W3 SMS 2 — Here we detect 27 molecular lines from 
11 species ( 12 CO, DCN, H 2 CO, HC 3 N, CH 3 OH, OCS, 
HNCO, SO, S0 2 , 13 CS, CH 3 CN) and 2 additional CO 
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Figu re 2. The SMA 1.3 mm continuum maps of W3 SMS1, SMS2 and SMS3 regions. Top panel: the VLA 22.5 GHz continuum image 
(Tieftrunk et al. 1997) overlaid with the SMA 1.3 mm continuum contours. The contour levels start at 5cr and are in steps of 5cr for SMS1 
(<7 =3.6 mjy beam -1 ), in steps of 4cr for SMS2 (a =2.5 mJy beam -1 ) and in steps of la in SMS3 region (<j =2.5 mjy beam -1 ). Bottom 
panels: the SMA 1.3 mm continuum maps of each region. The contour levels are the same as in the top panel. The dotted contours are 
negative features due to the missing flux with the same contours levels as the positive ones in each panel. Beams are shown at the bottom 
left of each panel. Filled triangles mark the millimeter sources detected. 
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Figure 3. Left panel: The LBT-LUCI X s -band image overlaid 
with the SMA 1.3 mm continuu m map of t he W3 SMS2 region. The 
Ks-band image is adapted from Bik et al. (2011). Right panel: The 
IRAC 4.5 /im short exposure time image ( Ruch et al. 2007) overlaid 
with the SMA 1.3 mm continuum image. The IRAC post-bed data 
were processed with pipeline version S18.7.0 and downloaded from 
the Spitzer archive. The contour levels are in steps of 4cr starting at 
5cr (<j = 2.5 mjy beam -1 ). The dashed contours are the negative 
features due to the missing flux with the same contours levels as 
the positive ones in each panel. The beam of the SMA map is 
shown at the bottom left corner. 

isotopologues ( 13 CO and C 18 0) with lower energy lev- 
els Slower A between 5.3 K and 315 K (Table EJ). The 
hydrogen recombination line H30a is also detected to- 
wards SMS2-MM4, MM5 and MM6. Figure [3 shows the 
integrated line maps of all species (except the CO iso- 
topologues which will be discussed in Sect. I3.3j) . 

Figure [7] shows that the hot and dense gas tracers, such 
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Figure 4. Upper and lower sideband spectra extracted towards 
SMS1-MM1 and SMS3-MM1 with a resolution of 2 km s -1 per 
channel. The SO line and the 12 CO line in SMS1-MM1 are not fully 
plotted, which go to 84 K and 94 K, respectively. The small gap 
in each spectrum is an effect of the data reduction process, since 
each sideband is divided into two parts to be processed separately. 
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Figure 5. W3 SMS1 molecular line integrated intensity maps overlaid on the SMA 1.3 mm continuum emission in the background. 
Contour levels start at 3a and continue in steps of 2a. The a value for each transition is shown in the respective map in mjy beam -1 . 
The dotted contours are the negative features due to the missing flux with the same contour levels as the positive ones in each panel. The 
integrated velocity ranges are shown in the bottom left part of each panel in km s — 1 . The synthesized beams of the molecular line integrated 
intensity images are shown in the bottom right of each panel. The (0, 0) point in each panel is R.A. 02 h 25 m 40.68 s Dec. +62°05 / 51.5 // 
(J2000). 



as CH3CN and CH3OH, all peak on the SMA continuum 
source SMS2-MM2, which indicates the massive hot core 
nature of this continuum source. The H2CO line also 
peaks at SMS2-MM2 and exhibits extended emission to- 
wards SMS2-MM1, MM3 and MM7. Other detected lines 
at SMS2 all peak at SMS2-MM2 except for DCN and 
13 CS. The DCN line shows extended emission towards 
SMS2-MM1, MM3 and MM7, and the 13 CS line peaks at 
SMS2-MM3 and extends towards SMS2-MM1. All these 
features show that the continuum sources SMS2-MM1, 
MM3 and MM7 are at a younger evolutionary stage than 
SMS2-MM2. The detection of the H30a line emission 
towards SMS2-MM4, MM5 and MM6 shows again that 
these three continuum sources are part of the UCH11 re- 
gion W3 C. 

The HC3N and H2CO integrated intensity maps show 
an extra peak to the north of the SMA continuum source 
SMS2-MM1 (Fig. EJ), and one IRAC 8 /im point source 
is found there (Fig. [8]). We name this 8 jam point source 
as IRS4-C, which is not visible at shorter wavelengths 



and does not have a mm continuum counterpart. These 
features show the peculiar properties of IRS4-C. 

W3 SMS 3 — Far fewer lines were observed in the W3 
SMS3 region. We detect 7 lines from 5 species ( 12 CO, 
DCN, CH3OH, SO and 13 CS) and 2 additional CO iso- 
topologues ( 13 CO and C 18 0) with lower energy levels 
Slower between 5.3 K and 35 K (Table [6j). Compared 
to the other two regions, the excitation temperatures of 
the lines we detect in this region are much lower, which 
indicates that W3 SMS3 has a much lower temperature 
than the other two. Figure [9] shows the integrated line 
images of all species (except the t hree CO isotopologues 
which will be discussed in Sect. I3.3|) . The DCN line 
emission peaks at SMS3-MM2 and the 13 CS line peaks 
at SMS3-MM1 and MM2. The CH 3 OH and SO also peak 
towards the region between SMS3-MM1 and SMS3-MM2 
and exhibit extended emission in the direction of one of 
the continuum sources. Because the continuum emission 
is usually more compact than most molecular emission 
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Figure 6. W3 SMS1 sulfur-oxides line integrated intensity map overlaid on the SMA 1.3 mm continuum emission in the background. 
Contour levels start at 3cr and continue in steps of 3cr, except for the SO (65 — 54) panel, in which the contour step is 10cr. The a value for 
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from SO or CH3OH, missing flux is unlikely to explain 
the discrepancy. The extended gas emission in regions 
without detected continuum rather indicates a peculiar 
chemistry. The small number of lines detected in SMS3 
indicates that SMS3 is chemically younger than SMS1 
and SMS2. 

3.3. The molecular outflows 

Figure [10] shows the 12 CO channel map of the single 
dish data only. The vi sr of SMS1, SMS2 and SMS3 are 
—39.3 km s _1 , —42.8 km s _1 and —42.8 km s _1 , respec- 
tively. A red-shifted outflow component which extends 
towards the southwest is found associated with SMS1. 
However, the blue-shifted component is relatively con- 
fined around the continuum peak. We believe this out- 
flow is in the northeast-southwest (NE-SW) direction. In 
the SMS2 region, one blue-shifted emission structure to- 
wards the northwest and northeast which is associated 
with the SMA continuum source SMS2-MM2 was de- 
tected (from panel —52.8 km s _1 to —45.6 km s _1 in 
Fig. flQ]). However, the single-dish channel map does not 
show much red-shifted emission in this region. Similarly 
in the SMS3 region, the single-dish channel map does not 
show much high velocity emission. We combine the SMA 
data and 30 m data to study the outflow properties. 

W3 SMS1— The top panels in Fig. HI from left to 
right, show the 12 CO(2 — 1) outflow image of the W3 
SMS1 region derived from the SMA data only, com- 
bined SMA+30 m data and the IRAM 30 m data only. 



The bottom panel in Fig. [TT] sh ows the SiO(5— 4) out- 
flow map observed with SMA. Ilmai et all ((2000) and 
iRodon et al.l (|2008l ) observed multiple outflows which are 
associated with W3 IRS5, and Fig. [TT] also shows the 
complicated outflow environment. By comparing our 
outflow data with previous work which has been done 
for this region dClaussen et al. 1984; Mitch ell et aDl991l; 
Choi et al.l [T993: Hasegawa et al. 1994; Imai et al. 2000; 



Gibb et all 120071 : lR,odon et al.l 120081 ). we identify four 



outflows in this region and mark them with four lines 
in the SMA+30m panel of Fig. [TTJ Figure [12 shows the 
position- velocity (PV) diagrams of the 12 CO SMA+30m 
outflow observations. 

Outflow-a: The NE-SW outflow shown in Fig.[[0]is de- 
noted as outflow-a, which is also seen in all the 12 CO out- 
flow panels of Fig. [TT] The blue-shifted part of outflow- 
a is not as prominent as the red-shifted part. This 
outflow a nd its asymmetr i c stru c ture has also been de - 
tected bv iClaussen et aH (119841): [ Mitchel l'eTall (|1991h : 
iChoi et all (j!993f ) and lHasegawa et al.l ([19941 ) . Since the 
CO o utflow emission traces the entrained gas (|Arce et al.l 
120071 ). it can be significantly affected by the environ- 
ment. Therefore, the asymmetric structure of outflow-a 
could be due to different properties of the ambient gas. 
Outflow "Flow A" identifi ed with H2O maser observa- 
tions by Ilmai et all (|200Q[ ) is in a similar direction as 
outflow-a. We believe that the outflow "Flow A" and 
our outflow-a actually trace different parts of the same 
outflow. Outflow "Flow A" traces the inner small scale 
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Figure 7. W3 SMS2 molecular line integrated intensity images overlaid on the SMA 1.3 mm continuum emission in the background. 
Contour levels start at 3<t and continue in steps of 2cr except for the CH3CN and SO (65 — 54) panels, in which the contour step is 5cr/level 
and lOcr/level, respectively . The a value for each transition is shown in the respective map in mJy beam -1 . The dotted contours are the 
negative features due to the missing flux with the same contour levels as the positive ones in each panel. The integrated velocity ranges 
are shown in the bottom left part of each panel in km s — 1 . The synthesized beams of the molecular line integrated intensity images are 
shown in the bottom right of each panel. The (0, 0) point in each panel is R.A. 02 h 25 m 31.22 s Dec. +62°06 / 25.52 // (J2000). 
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Figure 8. The IRAC 8 /im short exposure time image (Ruch et al. 
l2007h overlaid with the SMA HC 3 N (left) and H 2 CO (right) inte- 
grated intensity images. The filled triangles mark the SMA mil- 
limeter sources we detect. The IRAC post-bed data were processed 
with pipeline version S18.7.0 and downloaded from the Spitzer 
archive. Contour levels start at 3a and continue in steps of 2cr. 
The a value for each transition is shown in the respective map in 
mJy beam -1 . The dotted contours are the negative features due 
to the missing flux with the same contour levels as the positive 
ones in each panel. The synthesized beams of the molecular line 
integrated intensity images are shown in the bottom right of each 
panel. 

structure and outflow-a traces the outer large scale struc- 
ture. The different directions between these two out- 
flows could be due to precession of this outflow caused 
by the relative motions of the driving source of out- 
flow "Flow A" with respect to that of "Flow B" found 



by llmai et all ((2000). However, this outflow does not 
show up in our SiO obse r vation s or the SiO observa- 
tions done by IRodon et al.l ([20081 ). The emission of these 
two molecules does not necessarily stem from exactly the 
same environment because unlike CO outflow emission, 
SiO emission traces s trong shocks in dense molecular gas 
([Schilke et al.l [199 7bf ) . The other three outflows can only 
be detected in the SMA and SMA+30m panel. The PV 
diagram of outflow-a (panel a, Fig. [12]) shows that the far 
end of the red-shifted part resembles the outflow Hubble- 
law with increasing velocity at a large r distance from the 
outflow center (e.g.. lArce 

Outflow-b: Outflow-b is more or less in the east- west 
(E-W) direction, and is also shown in the SiO out- 
flow map. The PV diagr am (panel b, Fi g. [T2]) shows 
a weak Hubble-law profile (|Arce et al.l 20071). This out- 



flow was first detected by SiO emission (Gibb et al. 2007) 
and is only detected in the better resolu tion SMA and 
SMA+30m maps. IRodon et all (|2008f ) claimed that 
one of their SiO outflows (their S iO-c) could repre - 
sent the SiO outflow detected by IGibb et al.l (|2007D 
(i.e. our outflow-b ). However, the outflow detected by 
IRodon et al.l (|2008[ ) has the red-shifted part in the east 
and the blue-shifted part in the west, which is exactly 
the opposite of our outflow-b. Thus, we believe this 
SiO outflow trace s a different outflow than the one in 
IRodon et all (|2008l ). 

Outflow-c: All the PV diagrams (Fig. [T2j) are domi- 
nated by one strong outflow which is aligned very close 
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to the line o f sigh t at or near the zero-offset point. 
IRodon et al.l (j2QQ8h also reported an outflow which is 
aligned very close to the line of sight. Along the line c 
in Fig. [TTJ the SiO outflow map shows an opposite red- 
blue direction compared to the CO outflow map in the 
vicinity of SMS1-MM1. Thus we suggest there is another 
outflow, outflow-c, which is aligned close to the line of 
sight. Line c in Fig. [TT] does not represent the direc- 
tion of outflow-c, we just choose it to produce the PV 
diagram (panel c in Fig. [T2]h which shows a symmet- 
ric velocity profile. This outfl ow is also de t ected in the 
SiO outflow map (Fig. E}. IRodon et all (j2QQ8h iden- 
tified two SiO outflows which are aligned close to the 
line of sight, and outflow-c is likely the joint contribu- 
tion of these two SiO outflows. Figure [6] shows that the 
emission from all the SO2 and its isotopologues peaks 
on the SMA continuum source SMS1-MM1. Since S0 2 
is kno wn to be enriched by shock interactions with out- 
flows (jBeuther et al.l [2009). the strong SO2 emission at 
all different energy levels could be due to the shock driven 
by outflow-c. 

Outflow-d: Outflow-d is roughly in the north-south di- 
rection, and is also shown in the SiO outflow map. This 
N-S outflow mat ches the water maser outflow found by 
llmai et all (j2000f ) and one of the outflows detected from 
prope r motions of compact radio sources bv l Wilson et al.l 
(2003|). 

The bottom panel of Fig. [TT] shows that the central 
part of the blue-shifted SiO outflow emission is offset 
to the SW direction by ~ 2" from the continuum peak 
of SMS1-MM1. This offs et could be produc ed by the 
outflow SiO-c identified bv lR,odon etaLl (|2QQ8h (the red- 



shifted part is in the E direction and the blue-shifted 
part is in t he W direction). T wo additional outflows 
identified bv lR,odon etHI (j2QQ8l ) (their SiO-d and SiO-e) 
are not found in our observations. In principle, we cannot 
exclude that this outflow is just a wide-angle outflow, 
however, from our point of view the diverse data and the 
observations mentioned above rather favour a multiple 
outflow scenario. 

The outflow map of SMS1 shows that the driving 
source of all the outflows we found here is the SMA con- 
tinuum source SMS1-MM1, w hich was reso l ved in to four 
1.4 mm continuum sources bv IRodon et al.l (j2008[ ). They 
claim the driving source of their outflow SiO-e may be 
the NH3 feature close to the projected center of SiO-e. 
This NH3 feature is associated with our SMA continuum 
source SMS1-MM4, and we do not see an outflow driven 
by this source. The existence of multiple outflows driven 
by SMS1-MM1 also confirms that SMS1-MM1 is a proto - 
Trapezium system as proposed by Megeath et al.l ((2005). 



W3 SMS2— The top panels in Fig. [13] show the W3 
SMS2 12 CO(2 - 1) outflow maps. The SMA and 
SMA+30m panels show that the blue-shifted outflow 
emission is strong and extended. Although the red- 
shifted outflow emission is relatively weak and confined 
to the SMA continuum source SMS2-MM2 region, it still 
has a peak intensity of 16<j (la = 30 mjy beam -1 ) in the 
SMA outflow map. In the 30 m only outflow map, the 
red-shifted outflow emission is smoothed out and only 
the strong blue-shifted emission remains. 

The bottom panel in Fig. [13] shows the SO (65 — 54) 
outflow map observed with the SMA. The red-shifted 
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Figure 10. The single dish 12 CO channel map with a spectral resolution of 2.4 km s — 1 . The contour levels start at 5a and continue 
in steps of 30cr (la = 0.2 K). The crosses mark the positions of the SMA continuum sources SMS1-MM1, SMS 2- MM 2, SMS3-MM1 and 
SMS3-MM2. The arrow in the bottom left panel marks the outflow-a we identify. The v\ sr of SMSI, SMS2 and SMS3 are -39.3 kins -1 , 
—42.8 km s _1 and —42.8 km s _1 . The beam of the single dish data is shown in the bottom left corner of each panel. 



SO emission coincides with the blue-shifted component 
SO-a and only peaks on SMS2-MM2, the same as the 
12 CO(2 - 1) red-shifted emission. The blue-shifted SO 
emission shows two lobes, one lobe peaks on the con- 
tinuum source SMS2-MM2 that we denote as SO-a, and 
the other one extends towards the northeast direction 
that we denote as SO-b. There is also blue-shifted 12 CO 
emission in the SO-b direction, but it is not as clear as 
SO-b. Since the UCHu region W3 C is next to SMS2- 
MM2, the ambient molecular gas of SMS2-MM2 could 
easily be blown away and affect the outflow by making 
the red- and blue-shifted outflow components asymmet- 
ric. We believe that SMS2-MM2 is the driving source of 
this outflow and the outflow is in the northeast-southwest 
direction. However, we cannot rule out the possibility of 
multiple outflows. 

The left panel of Fig. [14] shows the PV diagrams of the 
combined SMA+30m 12 CO(2-l) observations. The PV- 
plot cuts follow the straight lines shown in Fig. [13] top 
middle panel. The diagram is dominated by strong blue- 
shifted emission, and the red-shifted emission is quite 
weak. The right panel of Fig. [14] shows the PV diagram 
of the SMA SO(6 5 - 5 4 ) observations. The PV-plot cuts 
follow the straight lines shown in Fig. [13] bottom panel. 
The diagram shows that a pair of blue- and red-shifted 
emission features remain close to the center of the out- 
flow (i.e. the SMS2-MM2 peak position) while another 
blue-shifted lobe resembles the Hubble-law with increas- 



ing velocity at a larger distance from the outflow center 

(jArce et al.ll2007ft 

iRidge fc Moord (|2001l ) found a bipolar outflow associ- 
ated with W3 SMSI in the direction NW-SE. The red- 
shifted emission of this outflow peaks towards SMSI but 
the blue-shifted emission extends towards SMS2. We do 
not detect this outflow, however, we do detect outflows 
in SMS2 (Fig. [15]). Although the red-shifted emission is 
not as strong as the blue-shifted emission, it is prominent 
enough to show that the outflow is associated with SMS2- 
MM2. We believe tha t the blue-shifted outflow emission 
IRidge fc Moorel (|2001f ) detected is driven by SMS2-MM2 
not SMS1-MM1 (aka. W3 IRS5). Furthermore, there is 
a - 4 km s" 1 difference between SMSI and SMS2. One 
could easily mistake the CO emission at the ^i sr of SMS2 
as blue-shifted outflow emission. 

W3 SMS3— Fig. H51 shows the combined SMA and 30 m 
outflow images of W3 SMS3. All panels exhibit strong 
blue-shifted emission associated with SMS2. Some weak 
red- and blue-shifted features are detected toward SMS3, 
but we do not find any outflow signature associated with 
SMS3. Perhaps the outflow activity is weak at SMS3, 
and due to the interference from the strong emission at 
SMS2, we cannot recover the outflow emission. This 
again indicates that SMS3 is at an early evolutionary 
stage, with hardly any star formation activity. 

3.4. Rotational structures 



12 



Wang et al. 



W3 SMS 




- Blue [-92,-48] 

- Red [-29,25] 



510 Blue [-58,-44] 

Red [-34,-22] 




10 



5 
A R.A. 



-5 



-10 



Figure 11. W3 SMS1 outflow maps observed with the SMA and IRAM 30 m telescope. The grey scale in all panels is the SMA 1.3 mm 
continuum image. Top panels: The full and dashed contours show the red- and blue-shifted 12 CO (2—1) outflow emission, respectively. The 
top left panel presents the SMA data only map, the middle one is the combined SMA+30 m data and the right one is the IRAM 30 m data 
only map. The velocity-integration ranges are shown in the top right. All the CO emissi on c ontour levels start at 4cr and continue in steps 
of 8cr. The four solid lines in the SMA+30m panel show the PV-cuts presented in Fig. [12] Bottom panel: The full and dashed contours 
show the red- and blue-shifted SiO(5 — 4) outflow emission, respectively. All the SiO emission contour levels start at 4cr and continue in 
steps of 2cr. The a value for each outflow emission is shown in the respective map in mJy beam -1 , except for the 30 m data only map 
which has units of K km s — 1 . The dotted contours are the negative features due to the missing flux with the same contour levels as the 
positive ones in each panel. The synthesized beam of the outflow images is shown in the bottom left corner of each panel. The (0, 0) point 
in each panel is R.A. 02 h 25 m 40.68 s Dec. +62°05 / 51.5 // (J2000). 



We use S0 2 (22 T , 15 - 23 6 , 16 ), H 2 CO(3 , 3 - 2 0>2 ), H30a, 
CH 3 OH(4 2 ,2 - 3i, 2 )E and HC 3 N(24 - 23) to study the 
rotational properties of the three regions. We find rota- 
tion structures in SMS1 and SMS2, and do not find any 
rotational signatures in SMS3. 

W3 SMS1— In the velocity map of S0 2 (22 7 ,i 5 - 23 6 ,ie) 
(Fig. [T6]) . we see a clear velocity gradient in the 
northwest-southeast (NW-SE) direction which is perpen- 
dicular to the direction of outflow-a (Fig. [TT]) indicating 
that t he envelope o f SMS 1-MM1 may rotate in this direc- 
tion. IRodon et al.l (|2QQ8l ) also observed a similar veloc- 
ity gradient from a lower sulfur dioxide transition. The 
H 2 CO(3o,3 — 2o, 2 ) velocity map (Fig. fT6|) shows a larger 
structure with complicated velocity signatures. 

The H30a velocity map (Fig. [16]) shows a velocity gra- 
dient in the NE-SW direction and follows the distribution 
of the compact radio sources detected bv lTieftrunk et all 
(1997), which is nearly perpendicular to the direction of 
outflow-b. Hence the H30a emission could be tracing the 
rotation of the ionized gas, although we cannot exclude 
the possibility of an ionized jet. This velocity gradient is 
different from the one seen in S0 2 (227^5 — 236, i6)> and 
outflow-d is found in the similar direction but with an 
inverse red-blue signature (Fig. [TT]). 

The position-velocity diagrams of the 
S0 2 (22 7 ,i 5 - 23 6 ,i 6 ), H 2 CO(3 ,3 - 2 , 2 ) and H30a 



emission are shown in Fig. [TT] The corresponding 
PV cuts of each panel are shown in Fig. [16] The PV 
diagrams of S0 2 (22 7j i 5 - 23 6 ,i6) and H 2 CO(3 ,3 - 2o, 2 ) 
show that the rotational structure is not in Keplerian 
motion, hence it may be just a rotating and infalling 
core similar to the toroids or envelope described by 
iCesaroni et al.l ([2)07). The rotational envelope detected 
by our S0 2 (227 ? i5 — 236,i6) observation has a size of 
~ 10 4 AU at the given distance of 1.95 kpc. This 
indicates that the proto- Trapezium system shares one 
common rota t ing en velope. The velocity jump found by 
IRodon et all (l2QQ8h in S0 2 emission is also confirmed 
by our H 2 CO observations. 

The PV diagram of the H30a emission shows quite a 
small velocity gradient for the radio recombination line. 
We extract the spectrum towards the emission peak in 
Tabled] and derive the FWHM line width from the Gaus- 
sian fitting as ~ 9.2 km s _1 . For a typical HCH11 region 
with a temperature of 8 000 K, the line width due to 
thermal broadening would be 19.1 km s _1 . Hence, the 
H30a emission is likely tracing the gas associated with 
kinematic properties, e.g., rotation and/or outflow, but 
we cannot properly differentiate among the options. 

W3 SMS2— In the velocity map of CH 3 OH(4 2 , 2 - 3i )2 )E 
(Fig. [18]), we see a clear velocity gradient in the NE- 
SW direction, which is in a similar direction as the blue- 
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Figure 12. Position- velocity diagrams of W3 SMS1 for the 
12 CO SMA+30m outflow observations with a velocity resolution 
of 1.2 km s — 1 . The PV -cut of each panel is shown in the top- 
middle panel of Fig. [TT] as the same letter a, b, c and d. The 
contour levels are from 10 to 90% from the peak emission (30.3 Jy 
beam -1 for panel a, 26.6 Jy beam -1 for panel b, 28.3 Jy beam -1 
for panel c, 31.4 Jy beam -1 for panel d) in steps of 10%. The v\ ST 
at —39.3 km s — 1 and the central position are marked by vertical 
and horizontal dashed lines. The central position of panel a and b 
is the SMS1-MM1 source, and for panel c and d it is the red-shifted 
outflow peak at R.A. 02 h 25 m 40.61 s Dec. +62°05 / 51.7 // (J2000). 

shifted outflow SO-b (bottom panel, Fig. [13]). CH 3 OH 
is well known as a molecule that traces shocks, cores 
and masers in star forming regions ( Beuther et al. 2005a; 
Uargensen et al.ll2QQ4t iSobolev et al.ll2QQ7t ). therefore, we 
propose that this velocity gradient traces the near-side 
of outflow SO-b. Another velocity gradient that is in the 
NW-SE direction is shown in both the CH 3 OH(4 2 2 - 
3i j2 )E and HC 3 N(24 - 23) velocity maps (Fig. |TSJ). 
This velocity gradient is also perpendicular to the direc- 
tion of outflow SO-b (bottom panel, Fig. H3J). CH 3 OH 
and HC 3 N are well k nown low-mass disk tracers (e.g. 
IGoldsmith et all 11999( 1. hence the NW-SE velocity gra- 
dient suggests a rotational structure. 

The PV diagrams of the CH 3 OH(4 2 , 2 - 3i, 2 )E and 
HC 3 N(24 - 23) emission are shown in Fig. [fij The 
PV-cut of CH 3 OH(4 2 , 2 - 3i, 2 )E goes through the SMS2- 
MM2 dust continuum peak and is marked with line a in 
Fig. [T8j The PV diagram shows that the blue-shifted 
CH 3 OH(4 2 , 2 - 3i, 2 )E resembles the outflow Hubble- 
law with increasing velocity at a la rger distance from 
the outflow center (jArce et al.ll2007f ). The PV-cuts of 
CH 3 OH(4 2 , 2 - 3i 2 )E panel b and HC 3 N(24-23) are 
shown in Fig. Hgf [line b in the CH 3 OH(4 2 , 2 - 3i, 2 )E 
panel and the line in HC 3 N(24— 23) panel] are perpen- 
dicular to line a (the direction of the outflow SO-b). 
These two PV diagrams show that the rotational struc- 
ture is not in Keplerian motion. With a size of ~6 300 
AU (d=1.95 kpc), this structu re is similar to th e ro- 
tational structure described in I Wang et~aT] (|2011f ) and 



iFallscheer et al.l (|2009[ ). i.e. large scale structure with a 
non-Keplerian velocity gradient. We do not find rota- 
tional structure for SMS2-MM1, and no rotational sig- 
nature is found in SMS3 either. 

3.5. Temperature determinations 

The dense gas (n> 10 5 cm -3 ) tracer CH 3 CN 
can be used as a temperature determinator (e.g. 
lLorenfc Mundvl Il984h . For W3 SMS1, four k = 0...3 
lines of the CH 3 CN(12/e — Ilk) ladder were detected and 
they do not peak on any of the SMA continuum sources, 
so a spectrum was extracted toward the CH 3 CN inte- 
grated emission peak with only the compact configura- 
tion data. For W3 SMS2, seven k = 0...6 lines of the 
CH 3 CN(12/ C — life) ladder were detected and and they all 
peak on the SMA continuum source SMS2-MM2, so we 
extract a spectrum toward the SMS2-MM2 continuum 
peak with only the compact configuration data. Both 
spectra are shown in Fig. [20j We assume Local Ther- 
modynamic Equilibrium (LTE) and optically thin emis- 
sion, and construct the rotational diagrams following th e 
method outlined in Appendix B of IZhang et al.l (|l998a). 
The level populations Nj^ were calculated from the 
Gaussian fitting results of the spectra shown in Fig. [20j 
The derived rotational diagrams are shown in Fig. [2TJ 
along with the linear fit of all the detected k components 
which gives the value of T rot - 116±46 K for SMS1 and 
~ 140 ± 30 K for SMS2-MM2. Studies of lWilner et ill 
(1994) revealed that T rot derived with this method agrees 
well with those obtained from the Large Velocity Gradi- 
ent (LVG) calculations. The average line-width of the 
CH 3 CN(12 /c - ll fe ) spectra towards SMS1 (top panel, 
Fig. |20|) is ~ 1.6 km s _1 . Considering our spectral reso- 
lution is 1.2 km s _1 , the linewidth is not fully resolved. 
For a temperature of ~116 K, the thermal line- width 
of CH 3 CN is ~ 0.36 km s _1 . For hot cores with tem- 
peratures o f ^100 K, the average CH 3 CN line- width i s 
>5 km s- 1 (jZhang et al.lll998all2Q07l : lWang et al.ll2011f ). 
Even for some IRDCs the lin e- width is still >2.8 km s _1 
(Beu ther fc Sridh aran 2007). All these features indicate 
that the gas from which the CH 3 CN emission originated 
is warm and at the same time shows a low level of tur- 
bulence. To our knowledge this has rarely (or not at all) 
been seen before. 

H 2 CO is a well know kinetic temperature determina- 
tor (e.g. Mangum & Wootten 1993; Jansen et al. 1994; 
iMiihle et all 120071 : iWatanabe fc Mitch ell 2008). Since it 
is a slightly asymmetric rotor, transitions between en- 
ergy levels with different K are only collisionally ex- 
cited. Therefore, comparing the level populations with 
different K components from the AJ = 1 transitions 
gives an estimate of the kinematic t emperature of the 
medium (M angum fc Woottenl I1993D . The spectra of 
H 2 CO(3o, 3 — 2o, 2 ) and H 2 CO(3 2)2 — 2 2) i) were extracted 
towards the continuum peaks of SMS1-MM1/2 and 
SMS2-MM1/2/3/7 with the compact configuration data. 
The spectra were processed with CLASS and are shown 
in Fig. [221 We assume LTE and optically thin emis- 
sion, and calculate the temperature Tkj n of these sources 
follow ing the method described by iMangum fc Woottenl 
(1993). The temperature results are shown in Table [5l 
The high temperature of SMS2-MM7 (98±33 K) could 
be due to the fact that it is heated by the nearby infrared 
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Figure 13. W3 SMS2 outflow images observed with the SMA and IRAM 30 m telescope. The grey scale in all panels is the SMA 1.3 mm 
continuum image. Top panels: The full and dashed contours show the red- and blue-shifted 1 CO(2— 1) outflow emission. The top left 
panel presents the SMA data only map, the middle one is the combined SMA+30 m data and the right one is the IRAM 30 m data only 
map. The velocity-integration regimes are shown in the top right. All the 12 CO emission contour levels start at 4cr, for the blue-shifted 
outflow the contour step is 8a and for the red-shifted outflow it is 4cr. Bottom panel: The full and dashed contours show the red- and 
blue-shifted SO (65 — 54) outflow emission. All the SO emission contour levels start at Aa and continue in steps of 2a. The solid lines in the 
SMA+30m panel and the SO panel show the PV-cuts presented in Fig. 1141 The a value for each outflow emission is shown in the respective 
map in mJy beam -1 , except for the 30 m data only map which has units of K km s -1 . The dotted contours are the negative features due 
to the missing flux with the same contour levels as the positive ones in each panel. The synthesized beam of the outflow images is shown 
in the bottom right corner of each panel. The (0, 0) point in each panel is R.A. 02 h 25 m 31.22 s Dec. +62°06 / 25.52 // (J2000). 
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Temperature estimation results for SMS1 and SMS2. 
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Figure 14. Position-velocity diagram of W3 SMS2 for the 
12 CO (2-1) SMA+30m data and the SO data observations with 
a velocity resolution of 1.2 km s -1 . The PV-cut is in the NE-SW 
direction with a PA of ~ 50° E of N (the cut is shown in Fig. [13)) . 
The contour levels are from 10 to 90% from the peak emission (16.9 
Jy beam -1 for the CO panel and 2.6 Jy beam -1 for the SO panel) 
in steps of 10%. The i;i sr at —42.8 km s -1 and the central position 
(i.e. the SMS2-MM2 peak position) are marked by vertical and 
horizontal dashed lines. 

source IRS4-a (Fig.©. 

Despite the possibility that the H2CO emission is not 
really optically thin and it might be affected by the miss- 
ing flux, the Tkin we obtain from the H2CO intensity ra- 
tio agrees quite well with the T rot we obtain from the 
CH 3 CN(12 fe - ll fe ) line fit for SMS1. The CH 3 CN and 
CH3OH integrated intensity maps (Fig. [5]) indicate that 



a T r ot is obtained from the rotational diagram of the CH3CN 
lines. 

b Tj^in is obtained from the H2CO line intensity ratios. 
c The spectrum we use to estimate this temperature is extracted 
from the CH3CN emission peak, which is close to SMS1-MM2, 
so we put this number here. 

SMS1-MM2 has a higher temperature than SMS1-MM1, 
and this is supported by the T^m we obtain from the 
H2CO intensity ratio. 

The LVG m odel calculations done by 

iMangum fc Woottenl (p~993) show that tempera- 
tures determined from the H2CO(3o,3 — 20,2) and 
H2CO(32,2 — 22,1) line intensity ratios have a large 
dispersion depending on the volume density where the 
line ratio is less than 3.5, which corresponds to a temper- 
ature of 70 K. Since we do not have the volume density 
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Figure 15. W3 SMS3 outflow images observed with the SMA and IRAM 30 m telescope. The grey scale in all panels is the SMA 1.3 mm 
continuum image. The full and dashed contours show the red- and blue-shifted 12 CO (2—1) outflow emission. The top left panel presents 
the SMA data only map, the middle one is the combined SMA+30 m data and the right one is the IRAM 30 m data only map. The 
velocity-integration regimes are shown in the top right. All the 12 CO emission contour levels start at 4cr, for the blue-shifted outflow the 
contour step is 8cr and for the red-shifted outflow it is Act. The a value for each outflow is shown in the respective map in mJy beam -1 , 
except for the 30 m data only map which has units of K km s — 1 . The dotted contours are the negative features due to the missing flux 
with the same contour levels as the positive ones in each panel. The synthesized beam of the outflow images is shown in the bottom right 
corner of each panel. The (0, 0) point in each panel is R.A. 02 h 25 m 29.49 s Dec. +62°06 / 00.6 // (J2000). 
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Figure 16. SO2 (227,15 — 236, 16), H2CO(3o,3 — 20,2) and H30ai velocity (1st) moment maps overlaid with the SMA 1.3 mm dust continuum 
of W3 SMS1. The contours start at 5cr and increase in steps of 5<r in all panels (lcr=3.6 mJy beam -1 ). The dotted contours are the 
negative features due to the missing flux wit h th e same contour levels as the positive ones in each panel. The dashed line s in each panel 
show the PV diagram cuts prese nted in Fig. [TT] The crosses and the letters mark the compact radio sources detected by Claussen et al. 
(1994) and Tieftrunk et al. (1997). All moment maps are clipped at the five sigma level of the respective line channel map. The synthesized 
beams of the moment maps are shown in the bottom left corner of each plot. The (0, 0) point in each panel is R.A. 02 h 25 m 40.68 s Dec. 
+62°05 / 51.5 // (J2000). 



information for our continuum sources, the temperature 
determination can only be done between 20 and 70 K 
for line ratios larger than 3.5. Therefore, the Tki n we 
obtain for SMS2-MM2 is not reliable. Nevertheless, the 
temperature we obtain from the CH3CN and H2CO 
measurements confirms the hot molecular core (HMC) 
nature of the continuum source SMS2-MM2. 

4. DISCUSSION 

4.1. Different evolutionary stages of the continuum 
sources 

The SMA and IRAM 30 m data together reveal three 
massive star formation regions in different evolution- 
ary stages. With our high resolution SMA observations 
(~2 900 AU at the given distance of 1.95 kpc) we find 
that all SCUBA continuum sources fragment into several 
cores. 

W3 SMS3 is considered to be the youngest re- 
gion since no radio o r infrared sources are found here 
(jMegeath et al.l 12008). Our results also support this, in 



particular the low number of molecular lines and the ab- 
sence of outflow emission. The large amount of con- 
tinuum flux filtered out in our observations indicates 
that the gas in this region is relatively uniformly dis- 
tributed. VLA observations of NH 3 (1,1) and NH 3 (2,2) 
found extended emission toward SMS3, and further- 
more, this extended emission was resolved into four NH3 
clumps surrounding t he core emission towards SMS3 
(jTieftrunk et al.lll998l ). Our continuum source SMS3- 
MM2 is associated with the main NH3 core, and MM1 
and MM3 are ass ociated with one of t he compact NH3 
clumps, clump 4 (jTieftrunk et al.lll998[ ). From the line- 
width of the single dish 13 CO(2 - 1 ) data (~4.9 km s " 1 ), 

we 



following the equation described in lPillai et al.l ()2C 
estimate the virial mass for SMS3 a s ~ 1563 M^, which is 
much larger than the virial mass Tieft runk et al.l (jl998) 
obtained (300-400 M ) from the NH3 observations and 
the mass we obtain from the SCUBA measurements (Ta- 
ble [TJ. Our 13 CO line may not be completely optically 
thin, and traces more or less all of the molecular gas, 
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Figure 17. Position- velocity diagrams derived for the cuts along the observed velocity gradient in Fig. [16] with a velocity resolution of 
1.2 km s _1 . The PV-cut for H30ai and SC>2(227 ? i5 — 236, 16) is marked in the corresponding 1st moment maps (Fig. [16]). The PV-cuts for 
H2CO(3o,3 — 2o,2) a and b are shown as line a and b in the H2CO(3o,3 — ^0,2) velocity map (Tig. 1160 . All the contour levels are from 10 to 
90% of the peak emission (264 mJy beam -1 for H30o;, 833 mJy beam -1 for H 2 CO(3 ,3 - 2 ,2) panel a, 525 mjy beam -1 H 2 CO(3 ,3 - 2 0j 2) 
panel b, 2.68 Jy beam -1 for SC>2(227 ? i5 — 236,16)) in steps of 10%. The offsets refer to the distance along the cuts from the dust continuum 



peak of SMS1-MM1. The v is 
horizontal dashed lines. 




at —39.3 km s and the central position (the SMS1-MM1 continuum peak) are marked by vertical and 



level sulfur-bearing molecular lines, peak on the strongest 
continuum source SMS1-MM1, which is associated with 
the HCH11 re g ion c luster W3 M and infrared source IRS 5. 
iRodon et al.l ([20Q8h resolved SMS1-MM1 into four 1.4 
mm continuum sources. Combining the strong outflows 
we find here, we can confirm that a massive multiple 
system is forming in SMS1-MM1. The fact that SMS1- 
MM2 is not associated with any NIR point sources and 
has far fewer lines than SMS1-MM1 is consistent with a 
starless core. The high temperature of SMS1-MM2 we 
estimate (Table [5]) might be due to the shock heating 
caused by the outflows driven by SMS1-MM1 (Fig. HTJl . 
SMS1-MM4 shows almost no lines aside from the three 
CO isotopologues , and it is also a s sociat ed with an NH3 
core detected by iTieftrunk et al.l (|1998h . Therefore we 
conclude that SMS1-MM4 might be at a very early evolu- 
tionary stage and a potential starless core. Since SMS1- 
MM4 is in the direction of the red-shifted component of 
outflow-b and there is SiO emission nearby, it might be 
a swept up ridge by the outflow. 



Figure 18. CH 3 OH(4 2 ,2 - 3i, 2 )E and HC 3 N(24 - 23) velocity 
(1st) moment maps overlaid with the SMA 1.3 mm dust continuum 
of W3 SMS2. The contours start at 5cr and increase in steps of 5cr 
in all panels (<j= 2.5 mJy beam -1 ). The dotted contours are the 
negative features due to the missing flux with the same contour 
levels as the positive ones in each panel. The dashed lines in each 
panel show the PV diagram cuts presented in Fig. [19] All moment 
maps were clipped at the 5cr level of the respective line channel 
map. The synthesized beams of the moment maps are shown in 
the bottom left corner of each plot. The (0, 0) point in each panel 
is R.A. 02 h 25 m 31.22 s Dec. +62°06 / 25.52 // (J2000). 



including the more diffuse envelope, whereas the NH3 
and SCUBA traces the dense gas in the cente r. There- 
fore, the NH3 virial mass ITieftrunk et~a71 ([1998) obtained 
should better represent the mass from SCUBA. Hence, 
SMS3 is approximately in virial equilibrium. Since the 
peak column density (Table [1]) is also above the pro- 
posed threshold for high-ma ss star formation of 1 g cm -2 
(jKrumholz <ST McKee 2008), we infer that SMS3 is at a 
very early stage of massive star formation. 

For the more evolv ed regions SMS 1 and SMS2 , both 
harbor NIR clusters (Oi ha et al.l 120041 iBik etHI l2011h 
and HCH11 regions (jTieftrunk et al.l 11997). and both 
show strong signatures of active star forming activity. 
However, the continuum sub-sources in these two regions 
all show different evolutionary features. 

Evolutionary sequence in SMS1 — Our high resolution ob- 
servations show a more complicated picture. Most of the 
molecular lines, especially some high excitation energy 



Evolutionary sequence in SMS2— Tief trunk et al.l (|l998) 
found a filamentary structure traced in NH3 starting at 
the HCH11 region W3 Ca and extending towards the 
northeast. Our SMA observations not only recover part 
of the filament but also detect continuum emission to- 
wards the UCH11 region W3 C, which is the most evolved 
source in SMS2. Most lines detected in SMS2 peak on 
the SMA continuum source SMS2-MM 2, which is asso- 
ciate d with the HCH11 region W3 Ca (jTieftrunk et al.l 
1199/1 ). We also detect outflows tow ards SMS2-MM2. The 
i nfrare d source IRS4 detected by I Wynn- Williams et al.l 
(|1972[ ) consists of two infrared sources, which are IRS4-a 
and IRS4-b (Fig. [3j). IRS4-a was classified as an 08- 
B0.5V star, and is consid ered to be the exciting source 
of UCH11 region W3 C (|Bik et al.ll201lh . The other IR 
source IRS4-b is associated with the SMA continuum 
source SMS2-MM2 and the HCH11 region, and it also 
dominates the mid- and far-infrared emission. All these 
features show that SMS2-MM2 is in the hot core evolu- 
tionary stage. It is surprising that the strongest contin- 
uum source in this region SMS2-MM1 shows only a few 
molecular lines and is in a younger evolutionary stage as a 
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Position- velocity diagrams derived for the cuts along the observed velocity gradient in Fig. [18] with a velocity resolution of 
The PV-cuts for the left and middle panel are marked as line a and b in the CHsOH(42,2 — 3i2_)E velocity map, respectively 
(Fig. [T6]) . The PV-cut for the HCsN(24 — 23) is marked as the straight line in its 1st moment map (Fig. fl8l ). All the contour levels are 
from 10 to 90% of the peak emission (from left to right, 621, 644 and 370 mJy beam -1 , respectively) in steps of 10%. The offsets refers 



Figure 19. 

1.2 km s -1 



to the distance along the cuts from the dust continuum peak of SMS2-MM2. The v\ s 
SMS 2- MM 2 peak position) are marked by vertical and horizontal dashed lines. 



at —42.8 km s and the central position (i.e. the 





Figure 20. Observed CH3CN(12fc — life) spectrum toward the 
CH3CN integrated intensity peak in SMS1 (top panel) and con- 
tinuum peak SMS2-MM2 (bottom panel) with only the compact 
configuration data. The dotted line shows the best Gaussian fit 
obtained with CLASS for the lines. 

potential massive starless core. Last but not least, SMS2- 
MM3 and SMS 2- MM 7 exhibit fewer lines than SMS2- 
MM1 (Table EJ) and have strong DCN emission (Fig. [7]), 
which makes them the youngest sources in SMS2. All 
these continuum sources show clear evolutionary stages 
from the most evolved UCH11 region to potential starless 
cores within 30 000 AU. 

In summary, at large scales, the whole W3 Main region 
harbors the relatively evolved Hn regions W3 A, W3 H 
and W 3 D (Fig. [2]), and a young cluster is forming inside. 
Our observations reveal three massive star forming sub- 
regions in different evolutionary stages. W3 SMS3 is in 
the youngest evolutionary stage, while W3 SMS1 and W3 
SMS2 are more evolved. Furthermore, within the subre- 
gions, we find the SMA continuum sources in different 
evolutionary stages, from potential starless cores to typ- 
ical high-mass hot cores. The data clearly show that the 



age difference of massive star formation regions not only 
exists on clump scales, but also within the clump. 

4.2. Astrochemistry 

Chemistry can be used to determine the evolution- 
ary stages of massive star- forming regions. The contin- 
uum sources show different chemical evolutionary stages. 
SMS2-MM2 shows typical HMC chemical characteristics, 
with many hot/dense gas tracer emission, e.g. CH3OH 
and CH 3 CN lines (Fig.0). The 13 CS and DCN emission 
offset from SMS2-MM2 also shows its HMC nature, since 
in a high temperature environment (~100 K) CS quickly 
react s with OH forming SO and SO2 (Beu ther et al.l 
2009), and DCN can be destr oyed rapidly by reaction 
with atomic hydrogen (| Wright et al .11 19961 : lHatchell et al.l 
1998). For the younger sources, i.e. the potential starless 
cores, DCN and 13 CS lines are also useful to determine 
evolutionary stages. 

DCN is reported as a indicator of the evolution of the 
dense cores, as it is liberated from the evaporated grain 
mantle and could be destroyed rap idly in the hot gas 
by reaction with at omic hydrogen (| Wright et al.l 1996; 
lHatchell et aTl[l99cl . The ratio of DCN/HCN is often 
used to determine the evolutionary stage of high-mass 
cores. Since we do not have the HCN observations, we 
cannot get this ratio to constrain the exact evolution- 
ary stages of the sources, however, the detection/non- 
detection of DCN provides us with a clue about the 
evolution. We detect DCN emission towards SMS2- 
MM1/3/7 and SMS3-MM2, which indicates that these 
four continuum sources are in an early evolutionary 
stage. The non-detection of DCN towards the starless 
cores SMS1-MM2 and SMS1-MM4 could be due to the 
fact that the jets driven by SMS1-MM1 heated the am- 
bient gas and destroyed the DCN molecules. We do not 
detect DCN towards SMS3-MM1 and SMS3-MM3 either, 
but 13 CS emission is detected towards these two sources, 
which indicates they are in very early evolutionary stage 
too. 

SMS1-MM1 also shows HMC chemical characteristics, 
with many emission lines and all the SO, SO2 isotopo- 
logues peaking on the continuum source (Figs. [5] and [6]). 
However, of the typical dense/ warm gas tracers CH3OH 
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Figure 21. LTE rotational diagrams of the CH 3 CN(12 fc - ll k ) 
lines shown in Fig. [20] The solid line is the linear fit to all 
the k comp o nents calculated following the equation (B6) from 
IZhang et all (!1998al V 

and CH3CN, only very weak CH3OH emission and no 
CH3CN emission is detected towards the continuum peak 
(Fig. |5j). The offset peak of CH 3 OH could be explained 
by the shock produced by the outflow, but we cannot ex- 
plain why the CH3CN peak is offset from the continuum 
peak. 

Another puzzle is that CH3OH lines show extended 
emission towards the area between very young sources 
SMS3-MM1 and SMS3-MM2 (Fig. [9]). These are peculiar 
chemical properties of the SMS3 region. 

4.3. Sequential and potential triggered star formation 

NIR imaging and spectroscopy studies by iBik et all 
(|2Q11[ ) show an age spread of at least 2-3 Myr for the 
massive st ars within the H11 regions in W3 Main. X-ray 
studies by Feigel son fc Townslevl ([20081 ) revealed a low- 
mass young stella r cluster with an age of >0.5 Myr and 
lOjha et aiT (|2004) obtained an age of the NIR cluster of 
0.3—1 Myr form the if s -band luminosity function. In 
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Figure 22. H2CO spectra towards the SMA 1.3 mm continuum 
sources SMS1-MM1/2 and SMS2-MM1/2/3/7 with only the com- 
pact configuration data. The dotted lines show the Gaussian fit 
obtained with CLASS for each spectrum. 

contrast, the UCH11 regions such as W3 C and W3 F 
only have an age of ^10 5 yr ([ Wood & Churchwell 1989; 
Churchwelll 120021 : iMottram et"all 120111 : iDavies et al.l 



201 If ), and our SMA continuum sources SMS1-MM1 and 



SMS2-MM2 which are associated with HCH11 regions are 
even younger (few times 10 4 yr, IMottram et al.l 120111 : 
IDavies et al.ll201lh . Furthermore, the potential starless 
cores in SMS3 are still at the onset of massive star for- 
mation. The age differences reveal different stellar pop- 
ulations. The OB stars associated with the diffuse Hn 
regions and the low-mass cluster members are the first 
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generation stars, and the exciting sources of the UCHii 
regions and our SMA continuum sources are the second 
generation stars. The second generation star formation 
might be triggered by the first generation stars or per- 
haps th e star formation activities l ast over a prolonged 
period (Feig elson fc Town slev 2008f ) . 

To investigate the interactions between the Hn re- 
gions and the molecular cloud, we construct the single 
dish 13 CO(2 - 1) channel map (Fig. [23]). The chan- 
nel map shows that the molecular cloud forms circle-like 
structures with velocity gradients around the Hn regions 
W3 A, H and D. For W3 A, in the panel -40 km s -1 of 
Fig.[23j the molecular cloud forms a "ring" around W3 A 
with two main clumps located at the northeast and west 
of the Hn region' s edge. These two main clumps move 
to the south a bit in the panel —38.8 km s _1 , but the 
"ring" is still there. In the panel —37 km s _1 , the "ring" 
is no longer present, and the clumps move to the south 
and join into one. Finally, in the panel —34 km s _1 , the 
clump moves to the southwest. These structures indi- 
cate a velocity gradient of the molecular cloud around 
the Hn region W3 A from the north to the south. A 
similar velocity gradient is also found around the Hn re- 
gion W3 H, which starts in the panel -44.8 km s _1 and 
the direction is from the east to the southwest. For the 
Hn region W3 D, a "ring" structure is found in the panel 
—43.6 km s _1 of Fig.[23j and the clumps move south and 
east of W3 D in the next three channels. These veloc- 
ity gradients suggest that the Hn regions are expanding 
and interacting with the ambient molecular gas. Further- 
more, the YSOs and potential starless cores of the SMS1, 
SMS2, and SMS3 regions are located close to the edge of 
the Hn regions. These facts together indicate that the 
Hn regions may have even triggered the star formation 
in these sub-regions. 

Similar structures with molecular gas distributed 
as a ring structure have als o been found around 
many Hii regions or bubb l es jP eharven g et al.l 120091 : 
iBeaumont fc Williams! 1 2 OlOl: [Fang et al. 2012|). 

In the S255 complex. I Wang et al.l ([2011[ ) revealed three 
massive star forming regions with different evolutionary 
stages within the same region. Based on the age differ- 
ence between the low-mass cluster and the central high- 
mass protostellar objects, they suggest that the low- mass 
stars form first under pressure of the Hn regions sitting 
around them, and then the low- to intermediate-mass 
stars may enhance the instability of the central high-mass 
cores and may potentially have triggered the formation 
of the central high-mass stars. Unlike the symmetric ge- 
ometry of the S255 complex, W3 Main is more randomly 
distributed. 

SMS1 sits between W3 A and W3 B, and is 
also associated with the dense young stellar cluster 
(jFeigelson fc Townslevll2008l : lMegeath et al .111 9961 ). which 
has the highest YSO density in W3 Main. The excit- 
ing s ource of W3 A has an age of ~3 Myr (jBik et al.l 
l2011f h and is much old er than the NIR cluster (0.3- 
1 Mvr. lOiha et al.l 120041 ). Therefore, W3 A may trigger 
the formation of the NIR cluster and even the central 
high - mass stars (IRS5). 

Tie ftrunk et al.l (jl998[ ) found a filament along the 
edge of W3 H and W3 D from NH3 observations. 
Tieft runk et al.l (jl998) suggested that this filament may 
be a swept-up layer of W3 H. As part of the filament, we 



suggest the SMA continuum sources in SMS2 may have 
formed under the influence of W3 D, and may even be 
triggered by W3 D in the dense gas at the head of the 
comet ary Hn region W3 D as discussed by iHoare et al.l 
(2007). Additionally, we see an "age gradient", from 
the southwest to the northeast, with the objects getting 
younger and younger (Fig. [2]). 

5. SUMMARY 

Combining the single dish IRAM 30 m observations 
and the high resolution SMA observations, we character- 
ize the different evolutionary stages of massive star for- 
mation regions within the W3 Main complex. W3 SMS1, 
SMS2, and SMS3 show varying dynamical and chemical 
properties, which indicates that they are in different evo- 
lutionary stages. Even within each subregion, massive 
cores with different evolutionary stages are found, e.g. 
in SMS2, the SMA continuum sources show evolutionary 
stages from the most evolved UCHii region to potential 
starless cores within 30 000 AU. 

We find outflows in SMS1 and SMS2. The absence of 
an outflow in SMS3 indicates it is at the onset of massive 
star formation. The multiple outflows in SMS1 confirm 
a multiple system is forming. 

We detect 36 molecular lines in SMS1, 27 lines in SMS2 
and 7 lines in SMS3. The temperature is determined 
from the CH3CN /c-ladder rotational diagram and H2CO 
line intensity ratios for continuum sources in SMS1 and 
SMS2. 

Rotational structures are found for SMS1-MM1 and 
SMS2-MM2, and in SMS1-MM1 we find a velocity gradi- 
ent in the hydrogen recombination line H30<x This may 
trace rotation and/or outflow, but we cannot distinguish 
between these scenarios. 

Evidence for potential interactions between the molec- 
ular cloud and the Hn regions is found in the 13 CO (2 — 1) 
channel map, which may indicate triggered star forma- 
tion. 
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Table 6 

Observed lines at all continuum peafc 



[GHz] 


Transition 


-^lower 

[K] 




SMS1 




SMS2 








SMS3 




MM1 


MM2 MM3 MM4 MM1 


MM2 


MM3 MM4 


MM5 


MM6 MM7 


MM1 


MM2 


MM3 


LSB 


217.105 


SiO(5 - 4) 


21 




V V 
















217.239 


DCN(3 - 2) 


10 
















v 




217.830 


33 SO(5 6 , 6 - 4 5 ,s) 


24 




















218.222 


H 2 CO(3 ' 3 - 2 2) 


10 




V V 




yj 












218.325 


HC 3 N(24 - 23) 


121 




















218.440 


CH 3 OH(4 2 ,2 -3i, 2 )E 


35 




V V 


V 








v 


V 


v 


218.476 


H 2 CO(3 2 ,2 - 2 2 ,i) 


58 


V 


V V 


V 


V 












218.760 


H 2 CO(3 2 ,i -2 2>0 ) 


58 


V 


V V 


V 














218.903 


OCS(18 - 17) 


89 


V 


V 


V 














219.276 


S0 2 (227 ? i5 — 236,16) 


394 


V 


V 
















219.355 


34 S0 2 (lli,ii - 10o,io) 


50 




















219.466 


S0 2 (v2=l)(22 2 20 -22i 21 ) 


983 




















219.560 


'C 18 0(2-l) 


5.3 




V V V 


v 


v 




v 








219.798 


HNCO(10 ,io-9o,9) 


48 


V 


V 


V 














219.949 


SO(6 5 -5 4 ) 


24 


V 


V V V 


V 


V 




V 


V 


V 


V 


220.079 


CH 3 OH(8o,8-7i, 6 )E 


85 




V 


V 














220.165 


S0 2 (v 2 =l)(16 3 ,i4 -16 2 ,i4) 


613 


V 


















220.399 


13 CO(2-l) 


5.3 


V 


V V V V 


V 


V 


V 


V V 


V 


V 


V 


220.585 


HNCO(10i, 9 -9i, 8 ) 


91 


V 




V 














220.594 


CH 3 CN(12 6 -11 6 ) 


315 






V 














220.620 


33 S0 2 (lli,n -10o,io) 


50 


V 


















220.641 


CH 3 CN(12 5 -11 5 ) 


237 




V 














220.679 


CH 3 CN(12 4 -11 4 ) 


173 






V 














220.709 


CH 3 CN(12 3 -11 3 ) 


123 




V V 


V 














220.730 


CH 3 CN(12 2 -11 2 ) 


87 




V 


V 














220.743 


CH 3 CN(12i-lli) 


65 




V 


V 














220.747 


CH 3 CN(12 -11 ) 


58 




V 


V 















USB 



229.348 


S0 2 (ll 5 ,7 -12 4>8 ) 


111 




V 


























229.545 


S0 2 (13 2 ,i2-13i,i3) 


828 


V 


























229.759 


CH 3 OH(8-i,8-7 ,7)E 


77 


V 






V 


V 


















229.858 


34 S0 2 (4 2 , 2 -3i,3) 


7.7 


V 


V 


























230.027 


CH 3 OH(3 2 , 2 -4 M )E 


28 


V 








V 


















230.538 


12 CO(2-l) 


5.5 


V 


V 




V 


V 




V 


V 


V 


V 


V 




V 


V 


230.965 


SO 2 (37i0, 2 8 — 389 j2 g) 


881 


V 




























231.061 


OCS(19-18) 


100 


V 










V 


















231.221 


13 CS(5o-4 ) 


22 




V 






V 




V 










V 


V 


V 


231.901 


H30a 




V 




V 










V 


V 


V 










231.981 


S0 2 (14 3 ,n - 142,12) 


593 


V 




























232.266 


S 18 0(5 6 - 4 5 ) 


37 


V 




























232.419 


CH 3 OH(10 2 ,8 -9 3>7 )A+ 


154 


V 


V 








V 


















232.946 


CH 3 OH(10_ 3 ,8 - ll_2,io)E 


179 












V 



















Note. — "^"indicates that this transition is detected towards this continuum source. "— " indicates that the spectrum extracted towards this source shows a pure 
negative feature due to the missing flux problem. 



